Unquenched large orbital magnetic moment in NiO 
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Magnetic properties of NiO are investigated by incorporating the spin-orbit interaction in the LSDA 
+ U scheme. It is found that the large part of orbital moment remains unquenched in NiO. The 
orbital moment contributes about fiL = 0.29/is to the total magnetic moment of M = 1.93/is, as 
leads to the orbital-to-spin angular momentum ratio of L/S = 0.36. The theoretical values are in 
good agreement with recent magnetic X-ray scattering measurements. 



PACS number: 71.20.Be, 71.70.Ej, 75.10.Lp, 75.50.Ee 



Electronic and magnetic structures of late 3c? transition 
metal (TM) mono-oxides (MnO, FeO, CoO, and NiO) 
have been extensively investigated over the last decades 
HI . Conventional band theory with the local spin density 
approximation (LSDA) fails to describe the electronic 
structures of the compounds. In the LSDA calculations, 
the energy gaps of MnO and NiO are underestimated . 
Even worse, FeO and CoO are predicted to be metallic 
which are, on the contrary, large gap insulators in nature. 
The main problem in the LSDA is the use of mean field- 
type exchange-correlation functional which is improper 
to describe localized 3d electrons. 

There have been several theoretical efforts to cure the 
deficiencies in the LSDA, for example, the self- interaction 
correction (SIC) scheme [§-|e|, the GW approximation 
(GWA) §,§, and the LSDA + U method §§. The 
SIC-LSDA is in the line of the extended LSDA by re- 
moving unphysical electron interaction with itself. In the 
GWA, the quasi-particle energy is obtained through the 
self-energy calculation to the lowest order in the screened 
Coulomb interaction W. The GWA method applied to 
NiO gives a rather good description of the energy gap 
size Jt|. Computational load, however, is very heavy in 
the GWA. The LSDA + U method overcomes the fail- 
ure of the LSDA by incorporating the on-site Coulomb 
correlation U of the multiband Hubbard model-like. In 
this method, localized 3d electrons are treated separately 
from delocalized sp electrons. As a result, all the 3d TM 
mono-oxides in the LSDA + U are obtained as insulators 
with well developed energy gaps which are comparable 
to experimental values ||. In this way, the energy gap 
problem in 3d TM mono-oxides is considered to be solved 
by various calculational methods. 

Due to the outermost characteristics of TM 3d elec- 
trons, atomic 3d orbitals are greatly deformed in solids 
by the crystal field and/or band hybridization effects. 
Hence, the orbital moment of 3d TM ion is usually 
quenched in solids, because it originates from atomic na- 
ture of involved atomic elements pQ ]. For example, all 
the 3d ferromagnetic transition metals of Fe, Co, and Ni 
show negligible orbital magnetic moments in the range of 
Mi ^ 0.1/xb 0]- -^ is, however, expected that the strong 
Coulomb correlation in 3d TM mono-oxides preserves the 
orbital moment of localized 3d electrons by reducing the 
ligand crystal field effects at metal ion sites. 



In fact, for CoO, it is easily conceived that the orbital 
moment is only partially quenched because the measured 
magnetic moment M — 3.4/j.b 0| simply exceeds the 
spin magnetic moment alone and the minority ti g band 
is occupied only two-thirds of its available states. There- 
fore, the existence of orbital moment in CoO has been 
stressed many times PJ^JlI . It is also shown that the 
LSDA + U method gives a good description of mag- 
netic structure of CoO with a large orbital moment of 
fXl ~ IjiB On the other hand, in the case of 

NiO, measured magnetic moments are in the range of 
M = 1.77 ~ 2.2 ± 0.2^ B |l|-0, which are comparable 
to the spin only moment of isolated Ni 2+ ion. Therefore, 
the magnetic moment in NiO has been fully attributed 
to the spin moment and the orbital moment is expected 
to be completely quenched. 

However, recent magnetic X-ray scattering measure- 
ment fTij] indicates that the orbital-to-spin angular mo- 
mentum ratio in NiO is as large as L/S = 0.34, far from 
fully quenched orbital moment. In the nonresonant mag- 
netic X-scattering, the separation of spin and orbital mo- 
ment is possible, because the spin and orbital moment 
densities have different geometrical prefactors in the scat- 
tering cross section that can be adjusted by changing ei- 
ther the scattering geometry or the X-ray polarization. In 
Ref. jl7| , the orbital moment in NiO is extracted by the 
polarization analysis of nonresonant magnetic-scattering 
intensities. This method has evidenced a large contribu- 
tion of the orbital moment to the total magnetic moment. 
They also found that the spin and orbital moments in 
NiO are collinear. 

Magnetic moment is one of the basic ground state 
quantities which should be provided by an appropri- 
ate band method. If the orbital moment is so large in 
NiO, all the previous attempts [p| JT^JT9[ | which tried to 
directly compare calculated spin moments with exper- 
imental moments are under a mistake. The linearized 
muffin-tin orbital (LMTO) calculation for NiO with an 
orbital-polarization (OP) correction (LSDA + OP) in a 
crystal-field basis yields the spin and orbital magnetic 
moments of (is = 1.43/x.b and (jlt, — 0.12^, respec- 
tively [^0|. Although the total magnetic moment is sig- 
nificantly improved to M — 1.55/xs with the OP correc- 
tion, as compared to M — 1.23/1 b from the LSDA, it 
is still smaller than the experimental measurements and 
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LATTICE CONSTANT (a.u.) 
FIG. 1. Magnetic moment behavior with varying the lat- 
tice constant in the LSDA. Inset is the behavior of the or- 
bital-to-spin angular momentum ratio L/S. Increase of all 
the magnetic moments is understood by the reduction of the 
crystal field and the band hybridization strengths with vol- 
ume expansion. 
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FIG. 2. Coulomb correlation effects on the magnetic mo- 
ments obtained by the LSDA + U method. We have used 
the exchange parameter of J = 0.89 eV. Inset is the or- 
bital-to-spin angular momentum ratio L/S. The total mag- 
netic moment M = 1.93/is and the ratio L/S = 0.36 for 
U — 8.0 eV are in good agreement with experimental values. 



the ratio of L/S = 0.17 is only half of the magnetic X-ray 
scattering value @. The SIC-LSDA method § yields 
a large orbital moment of [1l = 0.27/is, which, however, 
was suspected as a methodological artifact by other au- 
thors llqJLsf . In more recent study, the orbital moment is 
simply neglected JlSj ] . 

To determine the size of orbital moment in NiO, we 
have performed the LSDA + U calculations using the 
LMTO band method within the atomic sphere approxi- 
mation (ASA). The incorporation of the spin-orbit cou- 
pling into the LSDA + U method is known to give right 
orbital polarization in strongly correlated electron sys- 
tems Q . The LSDA + U Hamiltonian is given by 

^LSDA+(7 = "^LSDA — ^dc + "Hu (1) 

where the first term in the right hand side is the LSDA 
Hamiltonian and the second term is the double counting 
correction for the third term Hu- With the Coulomb 
interaction U and exchange interaction J parameters, one 
can write Tide and "Hu, respectively, as 

Hdc-if/A(A-l)-ij^A CT (A CT -l), (2) 

a 

H U = \Y1 V ( mm '' m " m "') n ^m" n ^ m "' 
{m},(7 

+ l [V(mm';m"m"')-V(mm';m"'m")} 

{m},cr 

* n mm" n m'm"' (3) 

where n G mml is the d occupation number matrix of spin 
cr and N a = Tr(n^ m ,), N = N+ + N~ . We relate the 
screened Coulomb interaction V(mm';m"m'") with the 
Slater integral F k ; 
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V{mm';m"m"') = ^c k (lm,lm")c k {lm' 1 lm"')F k , (4) 

fe=0 

where c k (Im, lm') is a Gaunt coefficient. For 3d electrons, 
three Slater integrals of F°,F 2 , and F 4 are involved in 
the calculation. Among those, the ratio of F 4 /F 2 is 
known to be constant ~ 0.625 for most 3d TM atoms 
pl| . Hence, the actual number of parameters is reduced 
from three F k, s to two of U and J, which are given by 
U = F° and J = (F 2 + F 4 )/14, respectively. To de- 
termine the orbital moment, the spin-orbit coupling is 
simultaneously included in the self-consistent variational 
loop |Q . We have assumed the antiferromagnetic order- 
ing state of type-II, and the experimental lattice constant 
of a = 7.893 a.u. is used. 

In Fig. |l|, magnetic moment behavior is shown with 
varying the lattice constant in the LSDA. Both the spin 
and orbital magnetic moments increase monotonically 
with increasing the lattice constant. This feature is un- 
derstandable in view of that the crystal field strength at 
Ni sites and the hybridization between Ni 3d and O 2p 
bands are reduced with volume expansion. At the exper- 
imental lattice constant of a = 7.893 a.u., the spin and 
orbital magnetic moments are obtained as = 1.12/i^ 
and fj,L = 0.15//B; respectively, which are consistent with 
existing results |20| . Both the spin and orbital mag- 
netic moments are larger in NiO than in fee Ni which 
has [is = 0.63fiB and = 0.06^s. This suggests that 
3d electrons are more localized in NiO than in Ni metal. 
Although the ratio of L/S = 0.27 in the LSDA is only 
slightly smaller than L/S = 0.34 in the experiment (see 
Table. the spin and orbital polarizations in the LSDA 
are not large enough because of underestimation of the 
Coulomb correlation between 3d electrons. 

Figure |^ shows the Coulomb correlation effects on the 
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TABLE I. Calculated spin (fis), orbital (hl), and total (M) magnetic moments in /j,b at the experimental lattice constant 
of a = 7.893 a.u.. L/S is the orbital-to-spin angular momentum ratio. 



Method 


LSDA a 


LSDA + OP 6 


SIC-LSDA C 


LSDA + U d 


Experiment 


MS 


1.12 


1.43 


1.53 


1.64 


1.90±0.2 e 


I'L 


0.15 


0.12 


0.27 


0.29 


0.32±0.05 e 


M 


1.27 


1.55 


1.80 


1.93 


1.77 / , 1.90 s , 2.2±0.2 e 


L/S 


0.27 


0.17 


0.35 


0.36 


0.34 e 



a Present results. 

b Reference |2Q|. 

c Reference ||JT 

d Present results with U 

e Reference Wfl. 

f Reference |JL5| . 

9 Reference 111 611. 



8.0 eV and J = 0.89 eV. 



magnetic moments obtained by the LSDA + U method. 
We have used the exchange parameter of J = 0.89 eV 
which is comparable to literature value ||. Once the 
strong Coulomb interaction is introduced between Ni 3d 
electrons, the magnetic moments increase substantially. 
The role of the Coulomb interaction is significant to lo- 
calize Ni 3d electrons. With U = 8.0 eV found in Ref. 
H and Ref. we have obtained the spin and orbital 
magnetic moments of fis = 1.64/ig and (1l — 0.29/ig, 
respectively. The total magnetic moment M = 1.93/1 s 
and the ratio L/S = 0.36 for U — 8.0 eV are in good 
agreement with the experimental data. The change of 
the Coulomb parameter U by 1.0 eV results in incre- 
ments of both the magnetic moments and the ratio L/S 
by A/15 ~ A/i L ~ 0.01/is and A(L/S) ~ 0.01, respec- 
tively. The overall total magnetic moment difference and 
the ratio L/S difference are only AM = 0.12/i# and 
A(L/S) = 0.05 inbetween U = 6.0 and 10 eV. Thus, 
all the magnetic moments and the ratio L/S are rather 
insensitive to the Coulomb interaction U within the em- 
ployed parameter range. 

Calculated magnetic moments are summarized in Ta- 
ble. | in comparison with experimental data and pre- 
vious results by other calculational methods. As men- 
tioned above, the sizes of magnetic moments are un- 
derestimated in the LSDA. The total magnetic moment 
of M = 1.27 fiB in the LSDA is much smaller than 
M = 1.77 ~ 2.2 ± 0.2/xg in experiments In the 

LSDA + OP, the spin moment is significantly improved 
to fig — 1-43/iB ^0|, which, however, is still smaller than 
the experimental values. More serious is the decrease of 
orbital moment in the LSDA + OP, which unfavorably 
makes the consistency between the theoretical and ex- 
perimental L/S ratio become worse. Both the SIC-LSDA 
(M = 1.80/i B ) H and the LSDA + U (M = 1.93/i B ) give 
values in agreement with experimental data. But it can 
be concluded that the LSDA + U results with U = 8.0 
eV and J = 0.89 eV, among various calculations, are 
especially in best agreement with experimental measure- 
ments. The ratio of L/S = 0.36 in the LSDA + U is 
also consistent with L/S = 0.34 in the magnetic X-ray 



scattering measurement 17 . 

In conclusion, we have found using the LSDA + U 
calculations that the orbital moment in NiO is not fully 
quenched, which is as large as \xl = 0.29/x_b. Both the 
total magnetic moment of M — 1.93/xb and the orbital- 
to-spin angular momentum ratio of L/S = 0.36 for U = 
8.0 eV are close to experimental values. The Coulomb 
correlation and the spin-orbit coupling are crucial to get 
right magnetic polarization in NiO. 
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